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ABSTRACT 



The CMB anisotropy on arc second range is examined to test the power 
spectrum of perturbations in the small scale region and, in particular, to 
estimate the mass of dominant dark matter particles. It is shown that for the 
simplest evolutionary history with standard recombination, three and four beam 
observations could discriminate the mass of dark matter particles in the interval 
O.bKeV < Mdm < ^KeV with an antenna beam (0.5 - 0.25) arc minute with 
amplitude 10~ 7 . This interval is the most interesting one for the problem of 
Large Scale Structure formation and evolution. 



Subject headings: cosmology: theory - cosmic microwave background - large 
scale structure of the Universe - dark matter 
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1. Introduction 

The investigation of the power spectrum of primordial perturbations is now a key 
problem of observational cosmology. The recent detection of large scale cosmic microwave 
background (CMB) anisotropics by COBE (Smoot et al. (1992)) has provided us with 
a unique opportunity to probe the properties of primordial perturbations in large and 
intermediate scales. These data allow to estimate the amplitude and slope of power spectra 
in the large scale area (Gorski et al. (1994), Bond (1994)). At the same time, the complex 
analysis of motion and spatial distribution of galaxies has provided valuable information 
about power spectrum at smaller scales (Peacock & Dodds (1994), da Costa et al. (1994)). 
These results show the large progress that has been made in observational cosmology during 
the last years. 

It can be expected that in the next years a strong effort will be concentrated around 
CMB anisotropics on scales from several arc-minutes to few degrees. Actually, these 
investigations could provide valuable information about the most fundamental cosmological 
parameters, namely the initial power spectrum, the Hubble constant, the cosmological 
constant, the baryon density, the ionisation history, the dynamic of recombination, etc. 

However, the huge potential of CMB investigations is not still exhausted with the 
problems mentioned above. In this paper we would like to proceed further and to reveal the 
close connection of small-scale CMB anisotropy with a fundamental open issue of modern 
cosmology that is the nature, mass and properties of the dominant dark matter particles 
or, in other words, the shape of the power spectrum on the smallest scales. This area of the 
spectrum is important in itself for a reconstruction of the evolution history of the universe, 
for a description of processes of galaxy formation and as an important test of inflation. 
And, what is perhaps especially interesting, the shape of the power spectrum in the short 
wave area is very sensitive to the mass of the dominant type of dark matter (DM) particles, 
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M DM - 

Particle physics gives us a set of possible candidates for DM particles. First of all 
there are three types of neutrinos, the expected masses of which are smaller than 10 eV. 
Then there is the massive long-lived majoron with mass in the KeV range (Akhmedov et 
al. (1993), Berezinsky & Valle (1993)). Finally, some fraction of DM density can be related 
to the neutralino with a mass which is more than 20 GeV. It is a very important task to 
discriminate among these versions and to single out the mass of the dominant type of DM 
particles. 

There are few possible ways to detect DM particles. The most direct method is, of 
course, the laboratory test of some of their possible properties (for example, with decays of 
particles). However, in this way we can only confirm (or reject) some expectations of the 
theory which indeed could be not so certain. A valid, though not simpler, alternative is 
the investigation of the power spectrum in the small scale region. This method provides a 
certain information about the mass of DM particles which is, of course, the most important 
parameter. The analysis of galaxies distribution cannot provide us with the necessary 
information because small scale perturbations disappeared during the nonlinear evolution 
of Large and Super Large Scale Structure (LSS and SLSS). However, some information 
about the small scale area of the power spectrum is retained in small-scale CMB anisotropy 
and can be, in principle, set off with three and four beam observations at angular scales of 
about 0.5 - 0.25 arc-minutes. 

The current more popular DM models are Cold Dark Matter (CDM, Blumenthal 
et al. (1984)) and Cold + Hot Dark Matter (CHDM, Klypin et al. (1993)). However, 
recent theoretical investigations of LSS formation and evolution (Doroshkevich (1995), 
Doroshkevich et al. (1995) hereafter DFGMM) show that the best description of these 
processes is provided by the Warm Dark Matter (WDM) model with the transfer 
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function proposed by Bardeen et al. (1986) (hereafter BBKS). It describes the continual 
transformation from Hot Dark Matter to CDM and allows to obtain typical scales of LSS 
and SLSS which are in best agreement with the observational estimates (Buryak et al. 
(1994), Doroshkevich & Tucker (1995)). The WDM model provides that the structures 
form during a reasonable redshift interval z ~ 5 — 20 for particles masses in the range 
Mdm ~ 0.5 — AKeV. Therefore the mass of DM particles predicted by the theory of 
structure formation is close to the possible mass of majoron, m ; , evaluated by Berezinsky & 
Valle (1993). However, the upper limit found in this paper, ra ; < 1.34:h 2 KeV (h = H o /100 
km/s/Mpc is a dimensionless Hubble constant), is smaller than the values expected from 
dynamical models. Although this problem needs to be investigated in more detail, a 
comprehensive analysis of some consequences of the WDM model becomes important. 

The differences between WDM and CDM models are significant only in the small scale 
area. On large and intermediate scales all the predictions of both models are the same. 
Therefore, the analysis of small-scale CMB anisotropy is a unique way to test the WDM 
model and to provide observational estimates of the mass of DM particles. 

In this paper we examine a very simple model which assumes a standard recombination 
at a redshift z rec = 1000 without secondary ionisation up to small redshifts. Our analysis 
is based on an approximate semi analytical description of the radiation transfer function 
because here our main aim is to test roughly the proposed approach. The numerical 
estimates listed below show that this way seems to be promising enough and deserves to 
be examined thoroughly. Further steps of the analysis can include on one hand a more 
comprehensive investigation of the process of LSS and SLSS formation which will allow to 
find a more profitable range of parameters of DM particles. On the other hand, detailed 
numerical simulations of CMB anisotropics could be performed for different cosmological 
scenarios. 
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The paper is organized as follows. In Sec. 2 we provide a short description of the theory 
of LSS and SLSS formation and evolution with WDM model in order to emphasize some of 
the principle differences between WDM and the scale-free CDM model. This description 
can be considered as a ground to prefer the first one as the most realistic DM model. In 
Sec. 3 small-scale CMB anisotropy for WDM and CDM models is investigated. In Sec. 4 
we discuss and summarize our conclusions. 

2. LSS and SLSS formation and evolution with a WDM model 

A quantitative, albeit approximate, description of the formation and evolution of LSS 
and SLSS has been recently proposed by Buryak et al. (1994), Doroshkevich (1995) and 
DFGMM. This theory starts from the Zel'dovich approximation (Zeldovich (1970)), the 
great potential of which has been repeatedly demonstrated, and shows that the evolution of 
LSS can be described by a few approximate expressions involving simply moments of the 
power spectrum. In particular, it establishes the relationship of characteristic measures of 
LSS and SLSS and of their evolutionary behavior with the power spectrum. Some aspects 
of the formation of SLSS have already been discussed by Demiahski & Doroshkevich (1992) 
and Buryak, Demiahski & Doroshkevich (1992) though the evolution of LSS could be more 
useful in discriminating the many DM models under study today. For instance, it allows to 
reject the first and simplest HDM model. 

Although there is little presently available information about LSS evolution, the 
acquisition of more and better quasar absorption spectra and their detailed and careful 
analysis provides an exciting possibility for the observation of the evolution with time of 
large-scale structure in the Universe up to redshifts z ~ 5. 

The problem of SLSS and LSS formation and evolution can be formulated in terms 
of typical scales or, in other words, in terms of moments of the power spectrum. BBKS 
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considered the set of even spectral moments, namely 

k 2 dk 



(k)k 2j (2.1) 



and used mainly three moments, <t which is identical to the amplitude of perturbations 
and the values 

< k 2 >= o\lol < k 4 >= o\lol (2.2) 

which play an important role in the process of structure formation. 

The expressions for <t and < k 2 > are divergent for a scale-invariant CDM spectrum. 
In reality, there is a cut-off at large k and according to WDM models (BBKS) this 
cut-off is directly related to the mass of the DM particles dominating the universe, i.e. 
k m ax Mdm{Mdm /Q-mh 2 ) 1 ^ 3 ■ For neutrino types of WDM models a similar mass 
dependence for the damping scale, k max oc Mdm, can be found (BBKS). 

In this paper we will concentrate our attention around three DM models with transfer 
functions 

Tcdm = /n(1 2 + 3 4 g 34g) (l + 3.89 9 + (16.1g) 2 + (5.46a) 3 + (6.71a) 4 ) (2.3a) 
Twdm = (1 + 1-7? + (4.3a) 3 / 2 + a 2 )" 1 x exp(-0.5oi? J - 0.5(oi? J ) 2 ) (2.36) 
T$ DM = (1 + 1-7? + (4.3o) 3 / 2 + q 2 )- 1 x exp(-0.5oi? m - 0.5(oi? m ) 2 ) (2.3c) 
q = k/(n m h 2 Mpc- 1 ) R f = 0.2(^f s i? m = 0.26^ 

where M^m is in KeV. The first model is simply the scale free CDM model, the second one 
is BBKS WDM model and the third one is a neutrino type WDM model. The two WDM 
models differ only for the mass dependence of the damping scale, namely Rj and R m . The 
numerical coefficients for the damping scales depend weakly on the particles composition 
and need to be corrected for various cosmological models. 
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Strictly speaking, the following theoretical analysis, based on the Zeldovich' theory, 
can be performed only for the WDM models and not for the scale free CDM model. This 
remark becomes clearly superfluous when referred to numerical simulations since in this 
case an artificial cut-off for the power spectrum is introduced. It is interesting that both 
theoretical predictions and simulated results for the DM spatial distribution are insensitive 
to the substitution of the CDM model (with a short wave cut-off) to WDM models. 
However, this is no longer valid for the spatial distribution of baryonic matter and galaxies 
(see, e.g., Fong et al. (1995)). 

Following Buryak et al. (1994) and DFGMM we will consider the distribution of 
LSS elements along a random straight line (core-sampling approach). In this manner 
the structure elements can be strictly defined both in theory (Zeldovich' pancakes) and 
in simulated and observational catalogues. The one dimensional spatial distribution of 
structure elements is found to be Poissonian. Therefore, the mean separation between the 
elements is the main physical characteristic of the structure and provides the quantitative 
description of structure evolution. Analysis performed by Buryak et al. (1994) and 
DFGMM shows that the mean ID separation of LSS elements rapidly decreases during the 
first period of the active structure creation. Later, the sequential merging of LSS elements 
leads to a slow increase of the mean separation. Here, we present without proof the main 
quantitative relations of the theory, referring the readers to the original papers for more 
details. 

According to Zeldovich' theory, the typical separation of structure elements along a 
random straight line, l an during the active phase of structure formation is defined through 
two first moments of the power spectrum as follows: 




y = (1 + z)/(l + z cr ) = (1 + z)V5/<7„, r 2 c = 3./ < k 2 > 



- 9 - 



where z is the redshift and z cr is the 'characteristic epoch' for the creation of structure 
elements (Doroshkevich (1984), Doroshkevich (1989)). The values z cr and r c are presented 
in Table I for CDM and WDM models. 

A more detailed theoretical analysis (DFGMM) shows that other even moments are 
equally important for structure evolution. Namely, for CDM like broad band power spectra 
with Harrison- Zeldovich asymptotic 

p(k)/27r 2 = (9/16)t ax kT 2 (k) (2.5) 

with a quadrupole anisotropy Q 2 = (17/2.73) X 10~ 6 and a horizon = GOOOh -1 Mpc } the 
LSS evolution due to random motion and merging can be approximately described by the 
following equation: 



1 

'1 + z) 2 H 2 Jo 



1di.™7TTZvT-2 I- P(k)G(kl dia )dk (2.6) 



G(x) = x I u exp( — u )du 
Jo 

where Z^'s again is a characteristic scale of the separation of LSS elements along a random 
straight line. A simpler approximation to this equation is 

V / o + ^(l + -)- 2 + ^' 

{'CO 

J- 2 = / kT 2 (k)dk 
Jo 

Thus, the structure evolution can be described through the typical scale l which is, of 
course, the -1 moment of the spectrum, 

<r_! = (9/16)C/ - 2 (2.7) 

Values of l for CDM and WDM models are listed in Table I. It is important that during 
the second evolutionary period (merging), the strong differences between the predictions of 
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CDM and WDM models are eliminated and the main quantitative results for the present 
matter distribution are the same (see Ilss(z = 0), Iq and L in Table I). 

For CDM and WDM^ 1 ), the LSS evolution according to relations (2.3) and (2.7) is 
shown in Fig.l where the values 



are presented vs. redshift z. The typical redshifts z m i n when l cr = l<u s are listed in Table I. 

A similar description can be given for the evolution of SLSS. In this case, the typical 
scale L is defined through the spectral moments <r_i and <t_ 2 (Buryak et al. (1994), 
DFGMM) and its values are presented in Table I for CDM and for the two versions of the 
WDM model. 

Numerous theoretical investigations and numerical simulations demonstrate clearly 
that the LSS is an intermediate asymptotic of the general process of formation and evolution 
of cosmic structures. Generally, structures are formed through rapid collapse of matter 
along one axis to form a Zel'dovich 'pancake' and with less rapid collapse along a second 
orthogonal axis to form a 'filament', that is, creating a 'ridge' of higher over-density within 
the pancake. They then evolve not only by merging together according to the (2.6), but also 
through a general motion of matter along the filaments towards the branch points of this 
'general network structure'. This latter process destroys the regular LSS elements, creating 
the sites for the formation of the richest clumps of matter. We refer the reader to Shandarin 
& Zeldovich (1989) for an excellent review of much of this work. Later these processes were 
discussed in more detail by Kofman et al. (1992) and by Buryak et al. (1994). 

Quantitative estimates of the structure evolution were tested recently with six 
numerical simulations performed for CDM and Broken Scale Invariance power spectra with 
various computational boxes (Doroshkevich et al. (1995)). These simulations cover an 
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interval 1 < (1 + z cr ) < 10 and, thus, allow to test various periods of structure evolution. 
This analysis shows that developing LSS can be found in the redshift interval 



At present, the farthest galaxies are found up to redshifts z ~ 3 — 4 (see, e.g., Giavalisco 
et al. (1995)) and regular LSS elements are observed in deep surveys (Buryak et al. (1994), 
Doroshkevich & Tucker (1995)). At the same time, the quiet evolution of LSS described by 
equation (2.6) is observed through the analysis of quasar absorption spectra (Doroshkevich 
& Turchaninov (1995)). These data allow to conclude that the quiet evolution of LSS took 
place at least for < z < 5 and from equation (2.9) we can infer the corresponding values 
of z cr . These are found to lie within the interval l + z cr sa2-^4 — 5 and correspond to 
WDM {1 ~> models with M DM « 0.5 - 2KeV for the transfer function (2.3b) and to WDM {2 ~> 
models with M^m ~ 0.5 — 3KeV for the transfer function (2.3c). 

A more detailed analysis shows some dispersion in the moments of structure creation 
that gives rise to some uncertainties in the estimates of z cr from (2.9) and M^m- For this 
reason, we investigate a wider range of values for Mdm, namely 0.1KeV< M^m < 8KeV. 
Additional theoretical investigations and corresponding numerical simulations need to be 
performed in order to narrow the intervals for z cr and M^m- 

A similar approach, based on the analysis of dynamical processes of LSS evolution, is 
an important part of theoretical cosmology allowing to employ the observational data for 
the spatial galaxy distribution together with numerical simulations. Through this method 
some important constraints can be set on possible parameters of DM particles and on the 
power spectrum. However, this is still an indirect way because any information about the 
actual processes of structure formation cannot be achieved even with a modern astronomical 
technique. At the present time, the farthest quasars can be found with a redshift z ~ 4.5 — 5 
only and even such information is a very poor one. Thus, to investigate these processes we 



0.25 < 



l + z 

1 + z a 



cr 



< 2.5 
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need to have some alternative ways which could provide a direct information about the 
small scale range of the power spectrum and, therefore, about the early processes of galaxy 
and structure formation. 

Such way is the investigation of CMB anisotropics on sub-minute scales. 

3. Small-scale CMB anisotropy 

The theoretical analysis of CMB anisotropics was fulfilled by many authors and now 
the general behavior of anisotropics is well known for most models. Detailed numerical 
computations have already been discussed in many papers (e.g. Peebles & Yu (1970), 
Wilson & Silk (1981), Bond & Efstathiou (1984), Bond & Efstathiou (1987), Vittorio & 
Silk (1984), White, Scott & Silk (1994)) and provide the best comparison of a specific model 
with observational data. However, as a first step towards new investigations of a model, 
an analytical (or semi analytical) method is sufficiently accurate and to be preferred. This 
kind of approach requires the use of an analytical fit for the radiation transfer function 
A ra d(k } fi) which relates the primordial power spectrum of matter perturbations to the 
CMB anisotropics. 

The simplest analytical expression for the radiation transfer function was proposed 
by Doroshkevich (1988) and Atrio-Barandela & Doroshkevich (1994). Here a modified 
version of the transfer function was employed. It takes into account the variation of metric 
parameters during the period of equality of relativistic and non relativistic matter density, 
as suggested by Hu & Sugiyama (1994), and provides a valid description of small and 
large-scale perturbations together with the first Doppler peak but decreases the amplitude of 
second and further peaks. Since our analysis is concentrated around small-scale anisotropy, 
this improved form of the transfer function is sufficiently accurate for our aims. A more 
accurate but significantly more cumbersome semi analytical approach has been proposed by 
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Hu & Sugiyama (1994). 

At present, the most popular method to describe CMB anisotropics is through the 
squares of the coefficients of spherical harmonics (the C/'s), which is very convenient for 
large and intermediate angular scales (see, e.g., Bond (1994)). However, for small-scale 
anisotropy the description with an angular correlation function is more suitable (for the 
subdegrees scale such approach was discussed in detail by Jorgensen et al. (1995)). It is 
interesting that it can also be performed in terms of spectral moments or typical scales. 
The specific damping of CMB anisotropics during the process of recombination requires to 
include in the analysis the odd moments as well. Thus, the coherence angle, 6 C} can be 
related to the cr_ 3 / 2 spectral moment and, what is more interesting for us, the expressions 
for a three and four beam experimental set up relate to cr_i/ 2 and <7i/ 2 respectively. 

3.1. The approximate description of small-scale CMB anisotropy 

As it is known (e.g. Bond & Efstathiou (1987), Efstathiou (1988)), the angular 
correlation function for the Harrison - Zeldovich primordial spectrum can be written as 
follows: 

C(9X) = f {^r 1 ) 2 ft T2 (k) f Q d^ rad e- s ^J (y) (3.1) 
y = 2kR h sin(9/2)(l - ji 2 ) 1/2 

s d = (kd rf i) 2 s s = (kR h e s f(i- fi 2 ) 

where 6 is the switch angle between the beams, T(k) is the usual matter transfer function, 
the radiation transfer function Al ad (k } //) describes the transformation of the power 
spectrum to the radiation field, Jo(x) is the Bessel function, the value d r 0.1R rec gives 
the damping of anisotropics during the recombination, and R rec are the present horizon 
and the horizon at the moment of recombination and we approximate the finite antenna 
beam size as a Gaussian of angular width 6 S . 
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The function A 2 rad (k, ji) can be written as 

A 2 rad (k,p) = A 2 (k) + 3 f i 2 B 2 (k) 
u 2 

A = — 1 -\ -(1 — Acosiu — u eq ) + 3cosu + 3a r ) (3.2a) 

6a r + u 

u 2 

B = -(6a r /u + Asiniu — u eq ) — 3sinu) (3.26) 

6a r + u 

with u = A;i?j|^ rec , a r = (4 / o 6 /3 / o 7 )(l + 2; rec ) _1 = 31/i 2 6 , u eq = kRj\ z=Zeq , 
z eq = Pml P-f = 40fi/i 2 2v ec where Ofc is the dimensionless baryon density, z is redshift 
and z rec = 1000 corresponds to the standard recombination. Clearly A 2 a(i 1 for small 
k and A 2 rad 9a 2 + 12.5(1 + 3/j 2 ) for large A;. The expression (3.2) matches the relation 
proposed by Atrio-Barandela & Doroshkevich (1994) for u eq = 0. Using hypergeometric 
functions we obtain a better agreement with numerical computations in the area of the first 
Doppler peak. For larger u the trigonometric approximation is sufficiently accurate. 

The function A^ ad (k } //) (3.2) does not include the diffusion damping of radiation which 
can be essential for small scale perturbations because the typical scale of such damping is 
dp 1.5(Ofc/i)~ 1/ ' 2 /i~ 1 Mpc. Let us note, however, that this factor can affect only some 
quantitative estimates but can not cancel the considered effect because here we have to deal 
with forced perturbations as the perturbations in the DM component generate perturbations 
in the baryons and CBR. Therefore, the diffusion of photons leads to a strong depression 
of the oscillatory terms in (3.2) and decreases the large k asymptotic of the function A ra d 
approximately up to 9a 2 and a strong dependence of the final results on appears. This 
effect is more important for small values of when the photon mean free-path increases 
strongly and the hydrodynamical treatment becomes incorrect in the smallest scales. In 
this case the kinetic analysis must be used and this effect needs to be considered further. 

The amplitude of the correlation function, C(0,6 S ) can be written as follows: 
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12 f Qrms\ f 00 

5 \ -to 
with S = Sd + S s . 



( \ ( — r f 2 fJ c \ ( d/j,Ar d exp( — S) — monopole — dipole (3-3) 
V T / Jo k Jo 



For small switch angles, 6 <C 6 C} the expression (3.1) can be rewritten in the following 
standard form: 

c(e,e s ) _ e 2 JL^_ JL\_ ( o 4) 

C(0,e s ) 2W 2 C 7l 4!#4 72 6!^ ^ • > 

where 6 C is the coherence angle and 71 and 72 are dimensionless factors describing the 

short-wave part of the power spectrum (see, e.g., Bond & Efstathiou (1987)). If damping 

effects are not taken into account and 6 S = 0, d r = in (3.1) then 71 is proportional to 

the zero spectral moment, <t , and is divergent for a CDM power spectrum. With the 

damping during recombination, 71 retains a finite value even for CDM power spectra and it 

is proportional to the moment (t_i/ 2 - A further expansion of the correlation function in the 

series for a standard CDM power spectrum can be performed only with a finite beam size 

(Bond & Efstathiou (1987)) since for 6 S — > 0, 72 — > 00. 

For WDM power spectra the exponential cut-off of the transfer function T(k) changes 
the factors 71 and 72 in (3.4) and the effect depends on the mass of DM particles, M^m- 
Therefore, this way can be used for the determination of the shape of the power spectrum 
in the short-wave area and, in principle, for the determination of M^m- For WDM spectra 
the factor 72 is convergent even for a point-like antenna. However, even in this case the 
beam smearing can have very important consequences. 



With two beam observations 

D 2 = < {Tl ~ T2)2 > = 2 (C(0, 9 a ) - C(9, 9.)) = C(0, e s )^ 2 (3.5a) 
1 o v c 

it is possible to find only the coherence angle 6 C and, thus, other experimental arrangements 
are more suitable for our aims. The three beam observations (Uson & Wilkinson (1984)) 
< (Tx + T 3 - 2T 2 f > 
r2 o 



^3 = A T2 ~ = 1.5C(0, 9 S ) - 2C(6 } S ) + 0.5C(2#, S ) 
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= c^e a )~ (3-56) 

are more sensitive to the small-scale region of the power spectrum and they allow to 
determine the factor 71. And, of course, the four beam observations 

D 4 = < (Ti + 3T 3 - 3T 2 - T 4 ) 2 > = 5C{ ^ ^ _ 7 5C ^ ^ + 3C , (2 ^ ^ _ Q 5C ^ 0s) 

= C(0,6 a )~ (3.5c) 
are the most effective way to test the possible mass of DM particles through the factor 72. 



3.2. The two and three beam observations 

For the further analysis it is convenient to introduce three characteristic lengths related 
to the transfer function as follows: 

l~} 3) = £ dkT 2 e~^ 2 (3.6a) 

/-_ 3 1} = dkk 2 T 2 e~^ 2 (3.66) 

CO 



/" 5 = / dkk 4 T 2 e~^r (3.6c) 



where d s = Rh6 s and we include in the lengths definition the smoothing effect of the 
antenna. For small values of 6 S (and d s ) the influence of this factor on /(_ 3 ) and is 
negligible and these lengths are well defined even for a point-like antenna. On the contrary 
the influence of the antenna is extremely important for the third length, /1, which is 
divergent for the scale-free CDM power spectrum and is very sensitive to the mass of DM 
particles, Mdm, for the WDM spectrum. 



The coherence angle can be written as (d s < d r for 6 S < 10' 

2 ±_ ^fd 2 ~ dj / ( _ 3) 

c ^ R 2 ° (90fM 2 ) 2 + 25/2 1 J 
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The amplitude of the correlation function Co and the characteristic length /(_ 3 ) are 



sensitive mainly to the comparatively large-scale range of the power spectrum where possible 
distortions of its primordial shape are not significant. Moreover, this interval of wave 
lengths can be described and investigated with more accuracy through the C/'s (see, e.g., 
Bond (1994)). This relation demonstrates the remarkable stability of the coherence angle 
relatively to various models because the Doppler peak range has a moderate influence on Co 
and /(_3). Therefore we can conclude that the coherence angle 6 C is mainly a function of the 
Hubble constant (as /(_ 3 ) oc h~ 2 ) and the damping factors so that two beam experiments 
promise to give important information about the damping during recombination and about 
possible processes of secondary ionisation. 

A three beam configuration seems to be more interesting for us as it is more sensitive 
to the small-scale range of power spectrum. Actually, 



The results of calculations of 71 for two WDM models together with 6 C are presented 
in Fig. 2 vs. the mass of DM particles, M^m for h = 0.5 and = 0.1. As it was expected, 
the coherence angle 6 C is a very weak function of M^m- At the same time, 71 increases of a 
factor three in the range of masses 0.1 — 8 KeV and this effect, in principle, can be set off 
with three beam observations. 

It is necessary to test the influence of the antenna damping as well. Our estimates 
of the antenna damping are also presented in Fig. 2 for three possible beam sizes, 
6 S = 0.5', 0.25' & 0.1'. They show that with an antenna angle of 6 S = 0.5' we are able to 
test only the upper limit of M^m < 0.5 KeV and more narrow beams need to be employed 
to investigate the most interesting range of 0.5 KeV < M^m < 2 — 4 KeV. The antenna 
angle 6 S = 0.25' seems to be sensitive enough up to M^m < 2 KeV whereas 6 S = 0.1' covers 




-3 




and with known coherence angle it allows to estimate the typical scale 
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the full interval of expected masses, M^m < 8 KeV . 



3.3. 



The four beam observations 



The four beam experiment had been proposed by Jorgensen et al. (1995) to investigate 



approach with a different geometrical configuration. 

Clearly a four beam experiment is further more sensitive to the small-scale range of the 
power spectrum and is further more suitable for the test of M^m- As it was noted above, 
for high order terms in the series (3.4) the coefficients become strong functions of the beam 
size. However, for the same beam size the four beam experiment is more sensitive to M^m 
as the factor 



include the characteristic length l\. 

In Fig. 3a and 3b we plot 72 as a function of the switch angle 6 for two beam sizes, 
6 S = 0.25' & 0.5' and for h = 0.5 and fi 6 = 0.1. We investigate two WDM models for 
different values of M^m within the interval 0.1 KeV < M^m < 8 KeV and the scale-free 
CDM model. 

These Figures demonstrate clearly the great potential of four beam observations for the 
study of small-scale regions of power spectra and, thus, for the detection of the mass of DM 
particles, M^m- As it could be expected, for small masses, M^m < 0.5 KeV, the function 
72 is approximately const, what is a natural result of the strong damping of power spectra 
in the short-wave region. However, as M^m grows the area of application of expressions 
(3.4) & (3.9) reduces. When M^m becomes > 1 KeV, the results of the computations for 



the sub degrees range of spectrum (the Doppler peak area). Here we propose to use such 




-5 



(3.9) 
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WDM models match the CDM case on a common asymptotic of the form, 

72 ~ ^0 c /9f (3.10) 

with 1 < n < 2. In this case expression (3.4) is no longer valid and 72 depends strongly on 
6. Nevertheless, the differences between the predictions of various models are large enough 
to be discriminated with narrow antenna beams. 

Fig. 3a and 3b demonstrate that the influence of the antenna beam is not very significant 
for 6 S 0.5', 0.25' and the upper limit of the mass range that can be investigated depends 
more on the model itself. For WDM' 1 ' the limit is set by M DM < 2 KeV while for WDM' 2 ' 
masses M^m ~ 4 KeV can be discriminated. Finally, it is possible to show that for both 
models the interval of masses M^m < 1 KeV can be investigated even with a beam size 

e s « v. 

These parameters clearly exceed the available ones. However, this is the most direct 
way to test the small-scale range of the power spectrum and, thus, to solve many key 
problems of modern cosmology. 

4. Discussion 

In this paper we attempt to demonstrate that the huge potential of CMB anisotropics 
is not exhausted with the most popular problems as the initial power spectrum, the Hubble 
constant, the cosmological constant, the baryon density, the ionisation history, the dynamic 
of recombination, etc. At present, sub-minute angular scales are not as popular as the 
sub-degree ones. However, our analysis shows the significant potential of investigating 
fluctuations at the smallest angular scales with three and especially four beam experiments. 

In the previous Sections we have stressed the close connection of small-scale CMB 
anisotropy with the nature and properties of dominant DM particles, and therefore with the 
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fundamental problem of the formation of various cosmic structures like galaxies, clusters of 
galaxies and elements of Large and Super Large Scale Structure. This way seems to be the 
most realistic to reveal many crucial problems of modern cosmology and modern particle 
physics because the nature of DM particles is one of the fundamental issues in which the 
interests of both branches of science are intersected. 

Our analysis has a rather preferential character. The main aim of this paper is to 
present under the simplest assumptions the principle abilities of the proposed approach 
rather than to give an exhaustive discussion of all the details. Thus, many important 
problems were only mentioned and undoubtedly require further investigations. The first of 
all is the effect of photon diffusion, which will change the quantitative estimates and can 
introduce a strong dependence of the results on the baryon density, fij. Special attention 
should also be devoted to the possibility of a non-standard recombination since secondary 
ionisation will strongly increase the damping scale. In these cases, the hydrodynamical 
approach used above is no longer valid (especially for small values of 0&) and the kinetic 
analysis needs to be performed. However, for intermediate values of (0& > 0.03) a more 
accurate analysis can only rectify but certainly not lose the strong dependence of the results 
of four beam observations on the mass of dark matter particles, M^m- 

Random noise is especially important for narrow beam observations and was discussed 
in detail during the eighties. Since we do not touch this very specific and complicated topic, 
we refer to the excellent reviews of Partridge (1988), Readhead et al. (1992) and White, 
Scott & Silk (1994). 

The problem of DM particles is very important in itself and it is sufficient ground to 
develop further the analysis of CMB anisotropy on sub-minute scales. However, here we 
would also like to emphasize the strong connection of this problem with many key issues of 
modern cosmology, like the fundamental process of galaxy formation and many enigmatic 
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results of quasar observations. Moreover, it relates to the old problem of the evolutionary 
history and spatial distribution of various components of matter, namely luminous, baryonic 
and dark matter. Any progress in these directions requires the collection of more and more 
detailed information about the small-scale range of the power spectrum and, thus, about 
the mass and other properties of DM particles. 

In Section 2 some arguments were given to show this close connection through 
both qualitative and methodological considerations. In particular, it is significant that 
fundamental properties of LSS and SLSS as well as the main parameters of small-scale 
CMB anisotropy can be described in a similar way using moments or characteristic lengths 
of the power spectrum. Thus, the initial perturbations of gravitational potential provide 
large-scale CMB anisotropy and the galaxies concentration in SLSS elements. On the other 
hand, small-scale CMB anisotropy and the process of galaxy formation are described by the 
short-wave area of the power spectrum. This approach emphasizes the common physical 
ground and complementary character of various observational phenomena in different 
branches of cosmology. At the same time it shows that only a complete analysis of CMB 
anisotropy in a wide range of angular scales could provide a reliable solution of all the above 
mentioned problems. 

This paper was supported in part by Danmarks Grundforskningsfond through 
its support for an establishment of Theoretical Astrophysics Center. A.D. wishes to 
acknowledge support from the Center of Cosmo-Particle Physics, Moscow. 
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15.3 
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Table 1: The M^m is in KeV and all the distances are presented in (Jlti 2 ) 1 Mpc. 
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FIGURE 2. The coherence angle 9 C fare min) and the factor 
y 1 vs. mass of DM particle for WDM 11 ' (dashed line) and WDM (2) 
(solid line) for h = 0.5 and = 0.1. Dot lines show the antenna 
damping for different antenna beam-sizes, 6 S . 
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FIGURE 3a. The factor y z vs. the switch angle 6 for the scale- 
free CDM model (dotted) and for WDM (1) models with different 
masses of DM particles for h = 0.5 and Q b =0.1. Vertical dot lines 
show the antenna beam — size. The dashed — dot — dot and long- 
dashed lines correspond to M DM =4 KeV and the M DM =8 KeV. 
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FIGURE 3b. The factor j z vs. the switch angle 8 for the scale- 
free CDM model (dotted) and for WDM (2) models with different 
masses of DM particles for h = 0.5 and Q b =0.1. Vertical dot lines 
show the antenna beam-size. The long-dashed line corresponds 
to M nu =8 KeV. 



